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Introduction 

A continuum model is developed for the analysis and prediction of the dynamic response of 
multiphase panels that are subjected to impulsive loading. The multiphase medium consists of any 
number of different elastic anisotropic constituents including the case of voids. The impulsive loading is 
applied on the external boundaries of panel in the form of time-dependent tractions and/or displacements. 
Acoustic loading can thus be modeled as a time-dependent pressure wave. The derived continuum theory 
results in a system of second-order equations in time. This system represents the governing elastodynamic 
equations for the discretized regions of the multiphase medium, the interfacial continuity conditions 
between these regions, and the externally applied boundary conditions. This system is solved by an 
explicit step-by-step procedure in time, and the resulting time-dependent displacement and stress response 
can be recorded at any point within the multiphase medium. In addition, parts of the discretized geometry 
may contain fluids. Thus, the air surrounding the multiphase medium can be simulated, with the acoustic 
loading applied to the air rather than directly on the multiphase panel. This is useful as it more 
realistically simulated the impact of a pressure wave on the surface of the panel, as wave reflection can be 
captured. Various cases of porous absorbers can be modeled by employing this theory. Examples include 
limp absorbers, where the solid ligaments are so thin and compliant that they can be assumed to have 
near-zero stiffness; rigid absorbers, where the ligaments are so stiff that they can be assumed to be 
motionless (with very high stiffness); and poroelastic absorbers, where both the ligaments and the air 
within the considered domain participate in the wave propagation. In the present investigation, the latter 
case is mainly considered. 

Theoretical model development is first presented. Then, results are given for the structural/acoustic 
performance of several foam core sandwich panels. Finally, heterogeneous panels, consisting of a stiff 
phase and a very compliant phase, are examined to determine the effect of the arrangement of the phases 
on the panel acoustic behavior. 


Full Three-Dimensional Theory 

In the present section a three-dimensional continuum theory is developed for the analysis of wave 
propagation in multiphase materials. To this end, consider a three-dimensional domain defined by 0 < vq < 
D,0<x 2 <H, and 0 < x 3 < L. This region is divided into N a , TVp and N y subcells in the X\, x 2 and x 3 
directions, respectively. Each subcell (aPy), whose dimensions are d a , and l y (a = l,..., N a , P = 1,..., 
7Vp, y = 1,. . . N y ), can be filled by a distinct anisotropic elastic material (solid or fluid) or is left empty thus 
forming a pore, see Figure 1. 
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The derivation of this theory starts from the elastodynamic equations of motion 


_(“Py) _ ,,(aPy)_y (apy) 

"’ Y dt 2 ' ’ 


i,j = 1,2,3 


( 1 ) 


where Y 1 , p <a(!y) and t denote the stresses, displacements, mass density of the material in subcell 

(a(3y) and time. The stresses are related to the strains c^ 7 1 by the anisotropic Hooke's law 

a (ap y ) =c (a Py ) e (aPy ) 5 i,j,k, l = 1,2,3 < 2 > 


where C^f /) are the components of the 4th-order stiffness tensor of the material in the subcell. 

The displacement field in the subcell (o.Py) is approximated by a second-order expansion in the local 
coordinates x{ a) , x^' 1 and x 3 y) , the origins of which are located at the center of the subcells, as follows: 
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where the time-dependent ^{ooco > 0 = 1 ,2,3), are the volume-average displacements in the subcell which 


together with the higher-order time-dependent terms ; ( l,m,n ) = 1,2; must be determined. Thus, for 


example, Jf4“Py) denotes the coefficient in the expansion of t/ ; - Y } the quadratic term associated with the 
x 1 -direction. The resulting strain components in subcell (a[3y) are given by 
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By averaging Equation (1) over the volume of subcell (a(3y), the following three relations are 
obtained 


(7) 


where /^“ooo)’ ^(ooo) an d ^“ooo) which were originally introduced by Aboudi (1986, 1987) are given 
in terms of the surface-average of the tractions evaluated at the surfaces y{ a) =±d a / 2 , y ^ 1 =±h^/2 and 
V 3 Y) = ±/ y /2 , respectively. Thus 
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where the surface-average of the tractions are given by 
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(13) 


Next, by consecutively multiplying Equation (1) by x/ a) , jtj anc ' x 3 (y) and averaging the results over the 
volume of the subcell (aPy) and integrating by parts, the following three set of three relations follow 
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In these equations, the quantities /{“fjd), an d ^wooi) can be expressed in terms of the surface- 

average of tractions as follows 
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In addition, ^,y“odd) denote the volume-average of the stress aJ“ Py 1 over the subcell namely, 


(19) 
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Let the vector S[qoo) be defined by the following six components of 5^“,^ : 
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In terms of the reduced matrix representation C pq ; p,q = 1,. . .,6; of the stiffness tensor C, ; «, it can be easily 
verified by employing the strain components expressions in Equation (6) that its components are given by 



Next, by consecutively multiplying Equation (1) by x^ a \ .ly 1 '" and x 3 2(y) and averaging the result 

over the volume of the subcell (aPy) and integrating by parts, the following three set of three relations 
follow 
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In these equations, are the volume -average first moments of the stresses which are given by 
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It can be easily shown that the relevant components of the vector S ,l + m + n ^ 0, that is defined by 
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where p = l,... ,6. 

By substituting Equations (7) in (23), (24) and (25) one obtains, respectively, that 
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Hence substitution of Equations (3 1) to (33) in (7) yields 
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Just like the surface-average of tractions: t t , t, and t, which were defined by Equations 
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In the following, these surface-average displacements will be related to the microvariables W ^) ; 

(Imn) = 0,1,2; in the displacements expansions. To this end, by substituting the displacements expansion 
(3) to (5) in (35) to (37), the following relations are obtained 
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Subtraction of the pair of two equations in (38), (39) and (40) yields, respectively, that 
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Similarly, adding the pair of two equations in (38), (39) and (40) yields, respectively, that 
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Substitution of Equations (44) to (46) in (34) provides the following expression for the second time 
derivative ^)-(ooo) m terms of the stress moments and the second time derivatives of the surface-average 
displacements: 
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Similarly, the expressions for ^[“oo) , ^/( 020 ) anc * ^[002 j* can be obtained by substituting Equation (47) 
in (44) to (46) yielding respectively 
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With the established values of given above in terms of the second-order time derivatives of the 

surface-average displacements and volume -average stress, moments one can proceed and express the 
values of , Y^f/ml) anc * ^vumn) i n these terms as follows. Equations (31) and (14), (32) and (15) 
and (33) and (16) provide, respectively, that 
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2/;p 2/ y 


(55) 


AT 


(«Py) __ 

3/(001) _ 


P ( “ PY) / y d 2 


12 dt 2 


(2) +(aPy) 

U: 


(2) — (aPY) 

U: 


, o(aPy) 
a 3/(000) 


(56) 


From Equations (8) to (10) and (17) to (19), the surface-average tractions can be obtained as 


OlifaPy) 

h 


+ d a j(g Py) 


+ 7 


(aPy) 


1 /( 000 ) ^ 1 1 /( 100 ) 


(57) 


(2) -(aPy) 


h 


p 

2 


(apy) 

2 /( 000 ) 


, )(«Py) 
~ l ~* y 2/(010) 


(58) 


(3) ± (a p y ) 

h 


±^K 


(apy) 

3/(000) 


+ K 


(aPy) 

3/(001) 


(59) 


Hence, it is possible to express these surface-average tractions in terms of the second-order time 
derivatives of the surface-average displacements and volume -average stresses and stress moments by 

substituting in Equations (57) to (59) the expression of > ^ 2 Hhm) anc * ^lUmn) which are given by 
Equations (51) to (56) yielding 


t/ = 


± 


P (aPT) 4 d 2 

160 dt 2 


(9 + (aPy) 

Uj + 

V 


(!) — (aPv) A 
«/ 


(2) ~*"(aPy) 

+ U} 


(2) (apy) (^) + (apy) 
u i - Uj 


(3) “(otpy) 

11/ 


(60) 


+ 


p (aPT) (/ a d 2 
12 dt 2 


(0+( a py) 

«/ 


0)-(opy) 

«/ 


± 


^ c(aPy) 

4 d a * mm) 


Msl 
4 hj 


c(«Py) 

°2/(010) 


3<7 a y(apy) (-,(apy) 

. ,2 J 3/(001) ^ °l/(000) 
4/y 
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( 2 ) — (aPy) 

t, = 


(2) — (aPy) 

t; 


p <0tPy) /7p 

d 2 

160 

dt 2 

p (op % 

d 2 

12 

dt 2 

p (apr) / y 

d 2 

160 

dt 2 

p (apr) / y 

d 2 \ 


0)+(aPy) (l) - (aPy) 


^(2)+(aPy) (2) (aPy) (2)'*'(aPy) (2) (aPy) 

Uj + Uj - Uj - Uj 


12 dr 


•) (l)-(aPy) (2)+(aPy) (2) — (aPy) /^(3) +(ccPy) (3)“( a py) 

- Uj - Uj - Uj + 7 Uj + Uj 

V 

_' 3l J aPt ' + _2Ec‘“Pv’ . 21 


Equations (60) to (62) can be written in the compact form 


•S) i(ooo) 


= [^](a PY ) d_ < 


_[ z ](a P y) 
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where a = 1,..., N a ; P = 1,..., Nf, y=l,... N y and A’"' 1 ’ 71 is a 18x18 matrix whose elements depend on the 
dimensions of the subcell (aPy) and the properties of the material filling this subcell, whereas L ( “ Py) is a 
matrix of the same dimensions whose elements depend on the dimensions of the subcell. 

The continuity of the interfacial tractions between neighboring subcells implies that 


0)+(aPy) 

(0 (a+l,p,y) 



(64) 

h 

= t t , i = 1,2,3, 

a = 1,. 

■-,N a -1, p = l,..W p , y = l ,...,N y 

(T "*~(aPy) 

(2) - (a,p,y+l) 



(65) 

h 

= ti , / = 1,2,3, 

a = 1,.. 

■,N a , P = l,...,A p , y = \,...,N y -1 

(3) "t~(aPy) 

A)“(a,p.y+1) 



(66) 

h 

= ti , / = 1,2,3, 

a = 1,.. 

,.,N a , P = l,...,Ap, y = \,...,N y -1 

In addition, the continuity of the displacements at the interfaces of the subcells yields 


0) +(«Py) 

(1) (a+l,p,y) 



(67) 

Ui 

, i 1, 2, 3, 

a = 1,. 

ii 

>- 

ii 

cn. 

7 

( 2 ) + (aPy) 

(2) — (a,p+l,y) 



(68) 

U i 

Wj- , i 1, 2, 3, 

a = 1,. 

■;N a , p = 1,..., Ap -1, y = l ,...,N y 

(2) + (aPy) 

(2) — (a,p,y+l) 



(69) 

Ui 

Wj , i 1, 2, 3, 

a = 1,.. 

;N a , P = l,...,Ap, y = \,...,N y — 1 

For imperfect bonding between the interfaces. Equations (67) to (69) take the form 



(O+taPr) (0— (a+l,p,T) 

0)+(aPy) 



u, - u, -R\f' 

(0 ti 

= 0, 

(70) 


i = 1,2,3, a = l 

,...,N a - 

■1, P = 1, • - -TV'p , y = l,...,A y 



(2) + (aPy) 

(2)-( a ,p+l,y) 

n(aPy) 

" A 2 i 

(2) +<aPy) 







Ui - 

- Ui - 

(0 t, 

= 

0, 




(71) 


i = 1,2,3, 

a = 1, 

..., A a , 

[3: 

=1,. 

,.,Ap 

-1, 

ii 


(3)+(aPy) 

( 3 ) _ (a,p,y+l) 

n(«Py) 

A 3( 

(2) +(aPy) 







“« - 

Ui - 

(0 ti 

= 

0, 




(72) 


i = 1,2,3, 

a = 1 

N 

a’ 

P 

=1, 

-Wp, 

Y 

= l,...,iV Y -l 



where (t ) , R^\t) and R^f >: ] (t) denote the time depending debonding functions that describe 

the behavior of the interfacial debonding of the subcell in the three directions. Wave propagation with 
imperfect bonding between the phases has been previously investigated by Aboudi (1988a). Not that, in 
the presence of perfect bonding, these functions are identically zero. 

Next, the following time-dependent boundary conditions must be imposed depending on whether the 
tractions or displacements are prescribed at the surfaces of the subcells at x x = 0, Xi= D; x 2 = 0, x 2 = H and 
x 3 = 0, x 3 = L. 
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and 


(D-dPy) 

t; = 




or 


(*> — (ipy) 

U: 




* = 1,2,3, P = l,..JV p , y - 


I, = 


0) + <iV o Py) / \ 

or u t = g- 1) (xf ) ,x 3 (r) ,t) i 

* = 1,2,3, p = l,..JV B , y = l,...,N y 


(73) 


(74) 


(2) (aly) 

t: = 




or 


* 2 ) (aly) 
**; 


i = 1,2,3, a = 


y = l,...,iV y 


(75) 


and 


Finally, 


and 


*2)+(oiVpy) / \ 


(2)+ ( ai 


(aiVpy) 


or u 
*=1,2,3, a = l,..jV a , y = 1,...,jV. 


= g F’k«‘ > ,*?>4 


Y 


(3) “(api) 

t: = 


/ \ (3) (api) / \ 

// 3> (* 1 ‘“ , ,*f,<) or = //’>(*<“>, jf,«) 

*' = 1,2,3, a = 1, - - -7V a , |3 = l,...,iVp 


(3)- 


’(apA*y ) 


= g/ 3) ( X l <a) , X 2 P) ,4 0r **,' =g,- 3) (4 a) ,^2 P) ,4 

* = 1,2,3, a = l,..JV a , P = l,...,iVp 


(76) 


(77) 


(78) 


In these equations, /- (i| and g- A) , i,k= 1,2,3, are time-dependent functions that describe the temporal 
form of the various applied loadings at the specific locations of boundaries of the region. 

2 0)±( a py) 2 ( 2 ) ± (a p T , 2 (3)± (a py) 

There are 1 8N a NnN y unknowns -4— «. , -4_ M , Ay M , z' = 1,2,3. On the other hand, the 

K dt 2 ' dr 1 dr 1 

interfacial traction and displacement provide 6(N a N$N y - NpN y ), 6(N a NpN y - N a N y ) and 6(N a NpN y - AyVp) 
which are respectively given by Equations (64) and (67) (or (70)), (65) and (68) (or (71)), and (66) and 
(69) (or (72)). In addition, the boundary conditions (73) and (74), (75) and (76) and (77) and (78) form 
another 6NpN y , 6N a N y and 6AyVp relations, respectively. Thus the total number of equations is 1 8N a NpN y . 
The resulting system of the 1 8N a NpN y equations can be formally represented by 

AX(0 = B(0 (79 


where A is a matrix of dimension 1 SN a NpN y x 1 8N a NpN y whose elements are the material parameters and 
geometric dimensions, and B(t) consists of time-dependent elements that involve the stresses and the 
imposed boundary conditions at time t. This equation represents the entire multiphase medium that result 
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from filling the subcells with various types of elastic materials (including pores and fluids). It also 
represents the elastodynamic equations that govern the motion of the material within the subcells, the 
various interfacial conditions between the subcells, and the applied time-dependent boundary conditions. 

A central finite difference of a second order accuracy in the time increment At reduces this ordinary 
differential equation to the following explicit form 

X(t + At) = 2 X(t) - X(t - At) + (At) 2 A _1 B(0 (80) 

from which the variables can be computed at time t + At from their known values at times t and t - At. 

This procedure is continued to the next time increment. Stability of this finite difference procedure is 
ensured by a proper choice of the value of the time increment At. It should be noted that in the analysis of 
porous materials, the tractions at the walls of the empty subcell are equal to zero, so that the number of 
equations can be reduced by excluding this type of subcell from the system of equations. 

In the special case of multi-layered composites, the accuracy of the model was verified by 
comparison with the ray theory prediction (Aboudi, 1988b). Furthermore, in this latter special situation of 
multi-layered composites, extensive applications and verifications of the accuracy of the model were 
recently presented by Clements et al. (1996, 1997a, 1997b, 1998) including comparisons with measured 
data. In addition, these authors introduced more developments and made more refinements to the model. 
Also, the analysis of harmonic wave propagation in laminated composites was presented by Aboudi and 
Hevroni (1991). 


Specialization to Two-Dimensional Theory for Thermoelastic Panels 

Consider the case of two-dimensional wave propagation such that there is no dependence of any field 
variable on the direction X\. This situation corresponds to the case in which the direction of the 
propagating waves is perpendicular to the Jti -direction. For example, wave propagation in fiber reinforced 
materials in which the waves propagate in the perpendicular direction to the continuous fibers can be 
analyzed by present two-dimensional theory. Wave propagation in a panel in the through-thickness 
direction can be analyzed by the present two-dimensional theory for wave traveling in the x 2 -direction. 

In this two-dimensional theory all variables depend only on the x 2 - and ^-directions, in addition to 
the time, t. The plane x 2 -x 3 is divided this time into Np x N y subcells with p = 1 , . . . , Np, y = 1 ,. . . N r In 
addition, x^ 1 and x]^ are local coordinates whose origin is located at the center of the subcell (Py), see 
Figure 2. Hence the displacement expansions, Equations (3) to (5), reduce in this special case to 


u 


(Py) _ , =(P WPy) , j(yWPy) 

' ■ A ") '' hint ' rr \ 


t 

+ — 
2 


1 ( 00 ) 

W )2 -^ 


V 


1 ( 10 ) 

/,2^ 

A 

4 


1 ( 01 ) 

f 


W, 


(Py) 


1 ( 20 ) 


+ — 


3x« 2 -^ 
v 3 4 v 


mAPy) 

" 1 ( 02 ) 


„<Py) _ Py) . APW(Py) + AyW(Py) 
"2 — " 2 ( 00 ) ' “*2 " 2 ( 10 ) " 2 ( 01 ) 


1 

H 

2 


f 


3x® 2 


v 


U 2 ^ 

_P 

4 


pfrPY) + J_ 
" 2 ( 20 ) + 9 


f 


1 2 


3x« 2 -fiL 

v 3 4y 


W. 


(Py) 


2 ( 02 ) 


,,(Py) _ rn(PY) i y(P)m(PY) + x ( y)u7(Py) 
“3 _ " 3 ( 00)^ x 2 " 3 ( 10 ) +a 3 " 3 ( 01 ) 


1 

+ — I 
2 


( 




3x‘P )2 -^ 




2\ 


3xP-^ 


m( Py) 

" 3 ( 02 ) 


V 


( 1 ) 


( 2 ) 


( 3 ) 
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The resulting strain components in subcell ( (3 7 ) are given by 


S<™ = 0 

F (Py) - m(Py) F ( P>tptPyt 
fc 22 ~ n 2(10)^ J 2 ''2(20) 

P (Py) - nHPy> + 3 r(yW(Py) 

fc 33 _ " 3 ( 01 ) ^ - u ’3 " 3 ( 02 ) 

g (py) _ J_(w(Py) + 3x (y) W (M + W^ y) + 3x^ ) W ^ y) ) 

b 23 2 V ^ 2 ( 01 ) ^ JA 3 " 2 ( 02 ) ^ " 3 ( 10 ) ^ JA 2 ^ 3 ( 20 ) / 

8 (py) _ jL(w(Py) + 33c (Y) wAPy) ) 
b 13 _ 2^ 1 (° 1 ) +JX3 " 1 ( 02 )/ 

F (py) _ l(ip(Py) , 3r ( P)fT ( Py) ) 
fc i2 - 2 V"i(io) + ja 2 "1(20)/ 

The governing equations that correspond to Equations (7), (15) and (16) are given by 

/<Pyi +/r(Py) _ n (Py)w(Py) 193 

^ 21(00) ^^31(00) _ K "l( 00)> l - 1,^,3 


AM 
“ 2 /( 10 ) 



* = 1,2,3 


tHPr) _ c(Py) — n (Py) -T_ii7(Py) 1 9 t. 

^ 3 /( 01 ) ° 3 i ( oo ) — P ,„"ir on ’ «— 


12 


<(01)’ 


Heterogeneous medium (panel) 



Subcell (py) 


Figure 2. — Geometry of the model specialized to two dimensions. 
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Here anc * ^mL,) are related to the surface-average tractions in the following manner 


Kpy) 

* / 2i(00) _ - 


(2)+ (flT ) (2) - (py) 

U - h 


( 8 ) 


jr(Py) _j_ 

^3/(00) - , 


and 




(py) 

2i(10) 


( 3 )+<Py) ( 3 )“(Py) 

h - h 


( 2 ) + (Py) ( 2 )“<Py) 

h + h 


(9) 


( 10 ) 


jvUh) -I 

^■31(01) _ 0 


where the surface-average tractions are given by 


( 3 )+(p T l t 3 ) — (Py) 

t : + t : 


( 2 )± 


(PO 1 f / T /2 


t , = — 


£ 


a 


(Py) 


/ y H/ 2 






(y) 


(ID 


( 12 ) 


(3) ^ = J_ f V 2 G (p/ 


-f 

Ap J-Ap/2 




v 




Furthermore, the volume -average stresses and stress moments are given by 

S t( mn) = = X ’ 2 ’ 3 


(13) 


(14) 


Let 


c(Py) _ c c c c c c "kpy) (15) 

mn ) \y\l(mn )’ ^22 (mn)’ *^33 (mn)’ ^23 (mn)’ *^13 (mn)’ ^\2{mn)\ 


In conjunction with strain components Equations (4), (14) yields 



2 r 


c(Py) -h_ 
>(10) - , 


/^(Py)jy(py) , /^(Py)rac(py) , WPyWfPy) 
L >6 >( 20 ) + >2 >( 20 ) + >4 >( 20 ) J 


(17) 
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7 2 

o(Py) _ Y 
^(01) - , 


c<? v« + c«} : + c®> V«g 


( 18 ) 


where p= 1,. . .,6 and T ( f y) are the components of the thermal stresses and 9 is the temperature deviation 
from a reference temperature, T R . 

The two equations that correspond to Equations (24) and (25) are presently given by 


hi 


hi 


_P_ 7-(Py) ,frP_ 
4 2 /( 00 ) + 12 


t'(Py) 
l iv 37(00) 


_ 9 c(Py) _ n (PY)^L 

za 27(10) _ P p 


ifr(PY) .^LmHPy) 

”7(00) + p ”7(20) 


i = 1,2, 3 


(19) 


/ 2 / 2 - 
IEHPy) ,2l 

4 ^37(00) + p 


r(PY) _ o e(Py) - n (PY)frL 
If 27(00) J Z,J 37(01) _ P 


/. 

12 


l 1 

mHPy) + 2 y_kHPy) 

"i(00) + 1Q ”7(02) 


i = 1,2,3 


( 20 ) 


Hence by substituting Equation (5) in these two equations, the following relations (that correspond to 
Equations (32) and (33)) can be presently established 


r(PY) _Hc(PY) , n (PY)^P_#(pY) 
J 2i(00) - ”27(10) + P 2Q ”7(20) 


( 21 ) 


r{ Py) _ c(Py) +n (PY)ilL#<PY) 
^37(00) - ”37(01) + P 20 "'(02) 


( 22 ) 


Hence, substituting Equations (21) and (22) in (5) yields 


1 ^_c(Py) . 1 ?.c(Py) -r,(PY) 
, 2 ”27(10) + ,2 ”37(01) “ P 

/7 P Y 




mHPy) __LmHPy) _ hEm(PY) 

”7(00) 2Q ”7(20) 2Q ”7(02) 


(23) 


(2) — (Py) ( 3 ) ± (Py) 

The surface-average displacements m, and w ; - are defined by 


( 2 ) ± (Pt) 

1 ( 

•ly/ 2 

«7 

=rJ 

i 

u 

-IJ2 ' 

(3) — (Py) 

i 

f*p/2 

«7 

^ ■ 

i: 

J -*p/2 


v(P) 


f 


: +- 


,(Py) 


x 3 w =±^ 
v 2 y 


dxi y) 




(24) 


(25) 


Substitution of the displacement expansions (1) to (3) in Equations (24) and (25) reveals that these 
surface-average displacements are related to the microvariables as follows 
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( 26 ) 


u - w^y) +_P_m7(Pt) . JLrodPt) 

u ; _ "i( 00) X 2 ^'(10) + ^ ";(20) 


(3)±(|)y) / 7 2 

M -ff(Py) +_Lo/(Py) + 22 m 7 (Py) 
u i "7(00) x 2 "((01) + ^ "((02) 


Manipulation of every pair in these equations results in the following 


. 1 + (Py) (“1 “(Py) 

W%1\=— Ul ~ U, 


wm=- 


1 R 3 )+(Pr) t 33 “(Py) 


Ui — lli 


, ? <2)+(lw (2) ( pt) 4 , 

»®i = T2 + », -Tl-W® 
"P J "P 


, 2 (3) (3) «*) 4 , 

u/(Pr>-_ u , ,, __wAPy) 

"((02) _ ,2 U i + U i ,2 "('( 00) 

T L J Y 

Substitution of Equations (30) and (3 1) in (23) yields 

1 d 2 |~ (2)+(|5r) (2)_ (^) (3)+ (im (3) “<py>l 60 1 m . l , R s 

mHPy) _ j, + ,, + M + ,, _i c'(Py) _i cCPY) 

"(•(00) - 14 dt 2 U i + M ( + U i + “(• + 7 p (Py) h 2 °27(10) + /2 °3((01) 

Consequently, Equations (28) to (31) establish the following expressions 


2 <7 2 r (2)+ <PO (2) “(Py>^| O+dSy) O-(Py) 240 1 , | 


2 <7 2 (2)+ *Py) (2 )- (Py , C (3)+ ( py) O) -(py)^ 240 ] , 1 , 

w no2) = ^j2~jp; ~ u i ~ u i + 6 Ui + »,• - 7 p (p y)/ 2 ^ S 2ilm + ~^ s m 

With the established values of , one obtains from Equations (21) and (6) that 

d n ( P y) d 2 T f (2)+ dW l 2 )-(Py)^ ( 3 )+(p r) (3)-(p y) l 79 19 ro 

° - a, - a, +2L^ ^Si S, 

' L v ) J 7/? p 11 ~y 


/-(Py) _ P ^ 
2,(00) 70 dt 2 
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( 2 ) + (Py) 

M, 


( 2 )“Py) 

Uj 


( 36 ) 


t-(Py) 

J 2/(10) 


P (P % d 2 

12 dt 2 


, c(Py) 

^ ° 2 /( 00 ) 


Similarly, Equations (22) and (7) result in 


K 


(Py) _ 

3/(00) _ 


P ( ^> d 2 
70 dt 2 


( 2 )+(Py) ( 2 ) — (Py) 

^( 3 ) +(py) 

Pl-dwV 

— Uj - Uj +6 

Uj 

+ Uj 


K 

) _ 


12 c(Py) , 72 c(Py) 

_,2 °2/(10) + _,2 3/(01) 

7//p 7/ y 


(37) 


and 


^(Py) _ P (Pt)/ y 
3i(01) ~ 12 <// 2 


(3)+(|, T ) (3) — ( py, 

«, - M; 


, c(Py) 

^ °3/(00) 


From Equations (8) to (11), the surface-average tractions can be obtained as 


( 2 > — (Py) 

h 


h 


+ P /(Py) , /ip 

x - ' y 2/(00) 2/(10) 


(Py) 


(38) 


(39) 


(3) — (Py) 

h 


±-*-k. 


(Py) 

3/(00) 


is( Py) 

' A 3/(01) 


(40) 


Consequently, Equations (35) to (38) yield that 


(2)±(Py) 

h 


= +- 


) (PY X d 2 


140 dt 2 

p cpY) /?p d 2 

12 ~dt 2 


^( 2 ) + (Py) 

Uj + 

V 

(2)+ (Pt , (2) — (py) 

U: ~ U ; 


( 2 ) “(py) 

U: 


\ < 3 > + (Py) 

(3) “(Py) 

1 

1 

Uj 

J 

- 

36 -(p Y ) 

6/; P C 

1\ 2im 

1 

"-4 1 
-^To| 
c 


3/(01) 


_!_ c(Py) 

^ °2i(00) 


(41) 


and 


(3) — (Py) 

P (Py V y 

d 2 

(2)+(Py, (2) — (Py) 

7 (3)+ ( py, (3)-,pyA 

tj 


- Uj - Uj + 6 

Uj + Uj 


140 

dt 2 


V )_ 


P (Pt> / 7 d 2 
12 dt 1 


(3)+,p y ) (3) “(py) 
U: ~ U, 


6/ Y_c(Py) | 36 e(p Y) 


Ik 


2 *-’2/(10) 
P 


11, 


■’3/(01) 


, c(Py) 

+ °3i(00) 


(42) 


These two sets of 12 equations provide the desired relations between the surface-average tractions and the 
second-order time derivatives of the surface-average displacements and the volume-average stresses and 
stress moments. 

These equations can summarized in the compact form 
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( 2 ) + 
w 3 

2 ( 2 )- 

- IV KPt) d u 3 

L lJ * 2 |< 3 > + 

Mj 

(3)- 


Pr) I J 33(00) 


where P = 1 iVp, y = 1 . . N y and [A' 1 ] l[iyl is a 12x12 matrix whose elements depend on the dimensions 
of the subcell (Py) and the properties of the material filling this subcell, whereas [Xi] <P/) is a matrix of the 
same dimensions whose elements depend on the dimensions of the subcell. In the present two- 

(2)±(Pt) (3)± ( |Jy) 

dimensional case there arc 12/V|>A, unknown surface-average displacements u t and u, 

The interfacial traction and displacements (assuming perfect bonding) between the subcells are 


(2)+ ( |3 T) 


(2) (P+l,y) 

: t; 


i = 1,2,3, 


P = l,...,Nn -1, y = 


(3) +(p T ) (3) — cp,y+i) 

t, = t : , 


1 = 1,2, 3, P = 1,...,7V b , y = l,...,iV -1 


(2)+ ( p y) (2) — CP+1, T ) 

U, = u { , i = 1,2,3, p = l,...,iVp-l, y = l,...,iV y 


(3) + (|5y) (3) -(p, r+1 ) 

Uj = Uj , / = 1,2,3, p = l,...,iVp, y = l,...,N y -1 


Finally, the boundary conditions at x 2 = 0 and x 2 = H are 


(2) (ly) / \ 


t 2 ) — (i T > / \ 

or u, 


i = 1,2,3, y = 
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and 


AtX3 


and 


( 2 ) + (A'py) 

h = 

*p»fe rt .() 

or 

( 2 ) + («py) 

Ui = 



i = 1,2,3, 

Y = 1, 



L the boundary conditions are 




(3>-(P„ 

h = 


or 

(3)-<pi) 

U, = j 



i = 1,2,3, 

P = l, 



< 3 > + (Pjv t ) 

h = 

4-”(4 w .4 

or 

<3) + ( pyv y ) 

1^ = 

sP’hf 


i = 1,2,3, 

P = l. 




( 49 ) 


(50) 


(51) 


In summary, just like the three-dimensional continuum theory, the interfacial tractions. Equations (44) 
and (45), and displacements, Equations (46) and (47), continuity conditions as well as the boundary 
conditions. Equations (48) to (51), provide the requested system of equations, given by Equation (79), of 
12AyV Y equations to be solved at any time yielding the surface-average displacements at that time. 


Extension to Inelastic Constituent Materials 

It is possible to analyze and model wave propagation in inelastic composites. Here it is assumed that 
the applied impact is sufficiently strong enough to cause the inelastic material to exhibit inelastic flow. To 
this end, the constitutive equation of the material filling subcell (aPy) is given instead of Equation (2) by 


t (“Py) _ £4c*Py) 


ijki 


p(aPy) 

e kl 


phaPr) 

~ S /W 


i,j,k, l = 1,2,3 


( 1 ) 


where s( ( a ^ 7 1 are the inelastic strain components. It is assumed that the inelasticity is governed by an 
isotropic flow rule. Thus for plasticity, the evolution equation of sjy “* !y) is given by 


gflaPr) _ ^(aPy) a (aPy) (2) 

where A (a * 3/) is the proportionality function. For viscoplastic constituents, the flow rule is given by 

gflaPy) _ yy(aPy) a (aPy) (3) 

Since the inelastic flow rule of the inelastic phase is isotropic, it must be that the stiffness tensor of 

the material is isotropic as well, namely 

c m = x««6 1( 5„ + nW’^s,, + a„a # ) < 4 > 

where /- <uflyl and p/' /flyl are the Lame' constants of the material. 

With the isotropic assumption of the inelastic flow rule. Equation (1) implies that 
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CT (apy) _ (^(aPY) s («Py) _2n( a Py)g / .( a Pl') 
ij ijkl kl ij 


( 5 ) 


In the perfectly elastic case, the quadratic displacement expansions. Equations (3) to (5), produce 
linear variation in strains and stresses at each point within a given subcell. In the presence of inelastic 
effects, however, a linear strain field generated by these expansions does not imply the linearity of the 
stress field due to the path-dependent deformation. Thus the displacement field microvariables must 
depend implicitly on the inelastic strain distributions, giving rise to a higher-order stress field than the 
linear strain field generated from the assumed displacement field representation. In the presence of 
inelastic effects, this higher order stress field is represented by a higher order Legendre polynomial 
expansion in the local coordinates. Therefore, the strain field generated from the assumed displacement 
field, and the resulting mechanical stress field, must be expressed in terms of Legendre polynomials as 

4“ w = 1 1 Z fc“’ k, fe w k fc”) 

1=0 m = 0 72=0 


00 00 00 / \ / \ / \ 

/ = 0 222 = 0 22 = 0 

where Y lmn ~ j(2l + \)(2m + \)(2n + I) , and the nondimensionalized variables Q’s defined in the interval 
-1 < Q < 1, are given in terms of the local subcell coordinates as C,\ a> = x (a) /(t/ a /2) , Cjf 1 = x (P y(/jp/2), 
/ 2 )- 

For the given displacement field representation the upper limits on the summation in the strains 
expansion (6) becomes 1. The upper limits on the summations in the stresses expansion (7) on the other 
hand are chosen so that an accurate representation of the stress fields is obtained within each subcell, 

which depends on the amount of inelastic flow. The coefficients anc * 1 !/(/»»!) i n the above 

expansions are determined as follows. 

The strain coefficients are explicitly determined in terms of the displacement field 

microvariables using orthogonal properties of Legendre polynomials. For example e { y“ood) ~ ^icmi anc * 

p(a(3y) _ JY A m(«Py) h 

*7/(100) ~ V J« a " 1 (2oo ) / z • 

The stress coefficients x ; y“^) are expressed in terms of strain coefficients and the u nk nown inelastic 

strain distributions by first substituting the Legendre polynomial representations for ej“^ y) and into 

the constitutive equations and then utilizing the orthogonality of Legendre polynomials as 

(aPy) _ ^(aPy) (aPy) _ n(aPy) (8) 

l ij(lmn) ^ ijko ^koilmn) Iv ij(lmn) 

The terms represent inelastic stress distributions and are calculated in the following manner: 


o(apy) _ .. 

A ij(Imn ) “ F((apY) 


r /„„, r 1 r 1 r 1 

4 J-iJ-iJ-i 


4 aWp i fe a) h fc P) h fe Y) Ki a, ^2 P) ^3 Y) 


(9) 
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With the above expression, the volume-average stresses and stress moments in the present case of 

inelastic constituents can be readily evaluated by employing Equation (26) in conjunction with (7) to (9). 
This yields the following expressions which respectively replace Equations (22) and (28) to (30): 


c(aPy) _ /-'’(apYW(apy) , /^(aPY)w(«PY) , r ( - a & ) W ( - a & ) 4- ) (m( a PY) , u/( a Py) ( 

°p(000) _ '-'pi "1(100) + ^p2 "2(010) "3(001) + ^ p4 2(001) + "3(010) / 

i -u w { - a & ) U (r^( a Py) + w (a Py) |_ n(aPy) 

+ ^p5 \"l(001) + "3(100) r ^p6 \"l(010) + "2(100)/ ^(000) 


( 10 ) 


■>(100) 


c(aPy) _ WapyW(aPy) WctPy)Ty(aPy) , Wa|3yW(apy) 

J n/iont — ^ ni moot T w "'tuonm 


p\ rr 1(200) ^ 6 rr 2(200) ^ 5 ^ 3(200) 


n(apy) 

2V3 " (100) 


(ID 


c(aPy) _ ^ 7 P [waPy W(«Py) , /^(aPy) w(apy) , ^(aPy)w(aPy) 1_ _^P_ ofaPy) 
">(010) - ^ [fr p6 "l(020) +l -p2 "2(020) ^ L y)4 "3(020)J A p(010) 


( 12 ) 


o(apy) _^[/^<apyW(apy) ^(aPyW(aPy) WaPyWtaPy) 1 ^Y p(aPy) 
">( 001 ) _ A lfr/>5 "l( 002 ) ^ ^4 "2(002>'->3 "3(002) J 2^3 ^t 001 


>(001) 


(13) 


where p = 1,. . .,6. The components of the vector are given by 


rj(aPy) _ [r, D 

**■( Imn ) l/M l(/m/2) ’ . 


p p p p ( a Py) 

22(lmn)> J ^33(lmn)’> JS ~23(lmn)’> ^13(/m«)? ^(/mrt). 


(14) 


These inelastic stress contributions are obtained by integrating the relevant flow rule to yield the inelastic 
strain at time t. 

In conclusion, the integration of the inelastic flow rule at the current time t provides the inelastic 
strain g /( “l 3y) from which the stress terms R j can be determined according to Equation (9). Hence, 

one can compute the volume -average stress and stress moments S from Equations (10) to (13). The 

1 9 — (aPy) ( 2 ) ± (aPr) ( 3 )± (a p T ) 

latter can be used in Equation (63) to obtain the various surface tractions t , t and t 
These can be employed in conjunction with the interfacial and boundary conditions to construct the 
current system of ordinary differential Equation (79) from which X at time t + At can be obtained. 


Structural/Acoustic Evaluation of Foam Sandwich Panels 

The above two-dimensional wave propagation theory has been applied to model the acoustic 
attenuation properties of structural panels of the type used in rotary aircraft cabins. These simulations 
were performed using a “virtual acoustic test chamber”, which defines the geometry in which the acoustic 
properties of the panels are evaluated. The main issue with a simulated test chamber is eliminating the 
effects of spurious reflections on the “measured” signal. For instance, in Figure 3, the applied cyclic 
loading, which simulates an acoustic signal, is applied on the lower boundary of the analysis space. The 
wave travels through the air and then impacts the bottom of the panel. The wave travels through the panel 
and into the air, where the wave is recorded at a “virtual sensor” location 0.5 cm above the panel surface. 
The wave then continues through the air to the boundary of the analysis space and is reflected back, 
eventually making its way back to the top of the panel. Since we are primarily interested in the effects of 
the panel on the propagating waves without reflected waves from the rear surfaces, we wish to exclude 
the effects of these latter waves. In order to avoid recording pressures that include the effects of these 
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reflected waves, the dimension of the virtual test chamber behind the panel must be sufficiently long. It 
was determined that a length of 21 cm was sufficiently long for an applied wave frequency of 10,000 Hz. 
Note that, as the frequency decreases, this length will increase as a longer time is needed for several 
waves to pass though the panel to be recorded, but this longer time requirement can allow reflected waves 
enough time to interfere. It should be noted that noise issues in many rotorcraft cabins tend to occur at 
frequencies in the range of 1000 to 3000 Hz. Examining these frequencies would require a longer virtual 
test chamber. Herein, air was assumed to be a fluid possessing no ability to transmit shear stresses whose 
bulk modulus is given by, 

K = pc/ (15) 

where c p is the velocity of sound in the fluid. The properties employed for air are, p = 0.00129 g/cm 3 and 
c p = 346 m/s. 

Figure 4 is a plot versus time of an applied 10,000 Hz signal (i.e., pressure wave) and the response at 
the location of the virtual pressure sensor (0.5 cm above the panel top surface location in the middle of the 
chamber) in the case of an empty chamber (filled only with air) and in the case of a typical panel. It takes 
approximately 250 ps for the wave to reach the virtual sensor location. Clearly, a large decrease in the 
wave amplitude is caused by the presence of the panel. 


Panel w 


Air 


Panel 


Air 


(a) 


TTTTT 

Applied cyclic loading 


21 cm 


Panel t 


7 cm 


CT21 =0 U2 - 0 



021 = 0 022 ~ applied 
023 - 0 


u-\ = 0 
U3 = 0 
O32 = 0 


(b) 


Figure 3. — Virtual acoustic test chamber used to simulate the acoustic behavior 
of structural panels using the wave propagation theory, (a) Dimensions. 

(b) Representative discretization and boundary conditions. 
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— Applied loading 
Signal in empty chamber 



Figure 4. — Example pressure waves recorded at the virtual 
sensor location 0.5 cm above the panel top surface. 

A facesheet thic kn ess study was performed on four different sandwich panels, with A1 and 1M7/8552 
facesheets and two Rohacell foam core densities, using the virtual test chamber (Figure 3). As shown in 
Figure 5, the total thickness of the panel was kept constant at 2 cm, while the facesheet (FS) thickness 
was varied between 0.8 and 1.6 mm. The relevant properties of the panel constituent materials are given 
in Table 1. Note that effective (homogenized) properties are used for the composite facesheet and foam 
materials, as opposed to modeling individual plies, fiber/matrix constituents, or pores within the foam. 
Figure 6 is a plot of sound pressure level (SPL) reduction caused by the presence of the panel in the 
virtual test chamber versus the panel areal mass for each of the four panel types with varying facesheet 
thickness. Note that, in this figure, the closer given panel design falls to the upper left, the better 
performance (as high SPL reduction with low mass is desirable). 


TABLE 1— MATERIAL PROPERTIES USED FOR 
PANEL ACOUSTIC ATTENUATION STUDIES 


Property 

Aluminum 

Quasi-Isotropic 

IM7/8552 

Rohacell 31A 
Foam 

Rohacell 110A 
Foam 

E, GPa 

72.4 

57 

0.036 

0.16 

V 

0.33 

0.32 

0.38 

0.38 

p, g/cm 3 

2.8 

1.578 

0.032 

0.11 


The SPL is given by, 


SPL = 20 log 


10 


f \ 
P rms 

Pref 


( 16 ) 


where p,ef= 20 pPa is the standard reference pressure for airborne sound, and the root mean square 
pressure is given by, 
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( 17 ) 


P rms 



av 

\[p(t)fdt 

0 


p{t ) is the pressure signal as a function of time, and t av is the duration of the pressure signal. 

Figure 6 shows that the denser Rohacell 110 foam core and the aluminum facesheet provide the best 
overall SPL reduction, but the composite facesheet with the Rohacell 110 provides almost as much SPL 
reduction at a significantly lower mass. For the Rohacell 110 foam core, increasing the FS thickness does 
not have much of effect on the SPL reduction. In contrast, the FS thic kn ess has a significant impact on the 
SPL reduction for the Rohacell 3 1 core as the contribution of this lower density foam to the panel 
damping is much lower than in the case of the Rohacell 110 foam core panel. It is noteworthy that the 
composite facesheet (FS) curves appear as almost an extension of the aluminum FS curves as if the 
aluminum facesheets continued to become thinner. 


2 cm 


1 

0.8 mm FS 




Figure 5. — Foam core panels 
considered with three different 
facesheet thicknesses. 
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Figure 6. — Areal mass versus SPL reduction for the four panel types simulated. 
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Figure 7. — (a) Panel bending/buckling performance index versus SPL reduction for the four panel types simulated, 
(b) Panel strength performance index versus SPL reduction for the four panel types simulated. 


From Figure 6, it appears that the Rohacell 110 core with thin composite facesheets provides the best 
overall performance of the panels considered. However, this figure considers only the raw unit mass of 
the panel and does not account for the panel’s structural performance. If instead we plot the SPL 
reduction versus structural performance indices (Ashby, 2005), as shown in Figure 7(a) and (b), it is clear 
that the less dense Rohacell 3 1 core panel with composite facesheets is also quite competitive. The 
performance indices are associated with the structural efficiency of the panel per unit mass for panel 
bending and buckling, Figure 7(a), and panel in-plane strength, Figure 7(b). The plots provide a good 
measure of the trade-off between acoustic performance and structural performance of the panels for 
lightweight applications. This is particularly important for rotorcraft cabin panels, where interior noise is 
an important design consideration as is structural mass. 

Acoustic Evaluation of Heterogeneous Panels 

Here we consider two panels consisting of 75 percent aluminum and 25 percent Rohacell 31A foam, 
as shown in Figure 8. The first includes two layers of Rohacell foam, while the second has discrete 
Rohacell 31 A inclusions. These panels were simulated in the virtual test chamber shown in Figure 3 at a 
variety of frequencies. The resulting SPL reductions are plotted versus frequency in Figure 9. Also plotted 
is the “mass law” for transmission loss through a panel (Fahy, 1985), 

TL = 201og 10 (in,/)-47 ( 18 ) 

where m s is the panel areal mass in kg/m 2 , /is the frequency in Hz, and TL is the transmission loss (i.e., 
SPL reduction) in dB. Both panels considered have an identical areal mass of 42.2 kg/m 2 ; all that is 
different is the arrangement of the foam and aluminum constituents. Equation (18) provides an 
engineering estimate of panel SPL reduction in frequencies governed by the mass law, where the material 
stiffness and microstructure do not play a roll. This is typically in the range of 1 00 to 1 000 Hz, as shown 
in Figure 10. 

Figure 9 shows that the panel with discrete foam inclusions provides more SPL reduction than does 
the layered panel. This is expected as the discrete inclusions are more effective in converting waves from 
longitudinal to shear. Both panels considered provide SPL reductions that are lower than the mass law 
prediction (as expected, see Figure 10). However, at the lower frequencies considered, the predictions 
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associated with both panel types converge and begin to follow the mass law curve. This is because, at 
lower frequencies, the waves have longer wavelengths and thus interact with the microstructure to a lesser 
extent. It can also be shown that the results are independent of the stiffness of the panel at these lower 
frequencies, as predicted by the mass law. 



0.5 cm 


Figure 8. — Heterogeneous Rohacell 31A/aluminum panels considered in the 
acoustic study. 



Figure 9. — Heterogeneous Rohacell 31 A/aluminum panel SPL reductions at a 
variety of frequencies. 
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Frequency, Hz 

Figure 10. — Typical frequency ranges for the acoustic attenuation properties of 
structural panels. 


Conclusion 

An integrated structural/acoustic model has been developed for the analysis of the dynamic response 
of multiphase panels that are subjected to impulsive loading, such as an acoustic pressure wave. Versions 
of the model applicable to three-dimensional bodies and two-dimensional bodies were presented. The 
geometry of the heterogeneous medium is discretized into an arbitrary number of subcells, each of which 
may contain a distinct solid material, fluid, or a void. By enforcing equilibrium and continuity conditions 
within and between these subcells, the dynamic response of the medium to time-dependent boundary 
conditions can be determined. An extension of the model to include inelastic material behavior was also 
discussed. 

The model was applied to analyze the structural/acoustic response of foam core sandwich panels, 
which are candidates for application in the cabins of rotorcraft. This was accomplished using a virtual test 
chamber, in which the acoustic wave can be applied to the air surrounding the panel rather than directly to 
the panel surface. This captures the interaction of the wave with the panel surface rather than assuming a 
known boundary condition at the panel surface, and is thus more physically representative of the 
interaction of the structure with the acoustic environment. However, as in real acoustic testing, care must 
be taken to ensure that reflected waves do not pollute the results so that the true panel acoustic 
performance is obtained. Various panel designs were compared in terms of their acoustic and structural 
performance, and it was demonstrated that a trade-off exists between the desire to design panels with high 
transmission loss (better acoustic performance) and the desire to design panels with the lowest mass 
(better structural performance). 

It is clear that the presented model is capable of evaluating the impact of material and geometric 
design parameter changes on the structural/acoustic performance of lightweight panels for application in 
rotorcraft. It’s ease of use and efficiency for performing full dynamic simulations, like those presented 
herein, make it an attractive alternative to the finite element method for this purpose. Enhancements that 
would make the model more useful would be an automated ability to alter the applied acoustic wave 
frequency, while simultaneously adjusting the size of the virtual test chamber to ensure that reflected 
waves do not influencing the results. 
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